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Abstract. Spectrophotometric observations of the coma of the comets 
C/2000 WM1 (LINEAR) and C/2002 Cl (Ikeya-Zhang) were taken during Nov, 
Dec 2001 and Mar, Apr 2002 respectively with 104-cm telescope of the State 
Observatory, Nainital. CN (A 3883 A) and C 2 swan bands ( A 4695, A 5165 
and A 5538 A) have been identified in both the comets. Na I emission was 
detected in comet Ikeya-Zhang. An estimate of CN and C 2 abundances and 
their production rates have been derived. Dust production rates have also been 
determined. 
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1. Introduction 

New comets come from the places in the solar system that are farthest and coldest and 
are believed to contain matter that is unchanged since the formation of the solar system. 
Therefore the study of constituents of the comets is important for understanding the early 
stages of solar system formation. At the center of the head of a comet is the nucleus. It is 
composed of chunks of matter. The most widely accepted theory of the comets, advanced 
in 1950 by Fred Whipple, is that the nucleus is ”a dirty snowball”. It constitutes ices of 
such molecules as water, carbon dioxide, ammonia, methane, with dust mixed in. This 
model explains many observed features of comets. The nucleus itself is too small to be 
observed directly from the earth. It is surrounded by coma which shines partly because 
of dust and gas reflecting sunlight. Gases liberated from the nucleus are excited by solar 
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ultra-violet radiation and radiate it in the visible, adding to the visibility of comets. The 
spectrum of a comet shows sets of lines from simple molecules. The comet is brightest 
and its tail is generally longest at the time of perihelion passage. However because of the 
orientation presented to the earth, it may not appear longest at this stage. As the comet 
recedes from the sun, its tail fades. 


Table 1. Basic data of the coma of the comets C/2000 WM1 (LINEAR) and C/2002 Cl 
(Ikcya-Zhang) at the time of observation 


Comet 

Date(UT) 

A 

(AU) 

r 

(AU) 

mi 

P 

(xl0 4 km) 

D 

(arcmin) 

C/2000 WM1 

Dec 

01.68, 2001 

0.317 

1.219 

5.9 

1.604 

2.33 

(LINEAR) 

Dec 

06.59, 2001 

0.329 

1.139 

5.7 

1.665 

2.33 

C/2002 Cl 

Mar 

23.57, 2002 

0.733 

0.519 

4.0 

3.588 

2.33 

(Ikeya-Zhang) 

Apr 

15.92, 2002 

0.451 

0.815 

4.9 

2.208 

2.33 


Apr 

16.92, 2002 

0.445 

0.830 

4.9 

2.188 

2.33 


Apr 

18.92, 2002 

0.434 

0.864 

5.1 

2.124 

2.33 


Apr 

29.90, 2002 

0.404 

1.048 

5.7 

1.977 

2.33 


Apr 

30.90, 2002 

0.405 

1.066 

5.8 

1.982 

2.33 

A = Geocentric distance; r = Heliocentric distance 

mi = Predicted integrated magntude; p = Radius of circular region in the sky at A 
D = Aperature diameter of the coma projected on the sky 


Table 2. Observed fluxes in 

emission bands of the comets for the observed dates of Table 1. 

Comet 

F(C 2 ,Av = 

0) 

F/F{ C 2 ,Av = 0) 



xlO 10 ergs / cm 2 / s CN(Av = 0) 

C 2 (Av = 1) 

C 2 (Av = 0) 

CN(Av =-l) 


5165 A 

3883A 

4695A 

5165A 

5538A 

C/2000 WM1 

1.165 

0.892 

0.109 

1.000 

0.114 

(LINEAR) 

2.346 

0.803 

0.024 

1.000 

0.247 

C/2002 Cl 

139.306 

1.292 

0.302 

1.000 

0.369 

(Ikeya-Zhang) 

62.486 

1.397 

0.008 

1.000 

0.379 


52.060 

1.265 

0.015 

1.000 

0.091 


27.598 

1.694 

0.146 

1.000 

0.615 


17.531 

1.756 

0.015 

1.000 

0.629 


00.548 

2.418 

0.058 

1.000 

1.000 


The comet C/2000 WM1 (LINEAR) was detected by LINEAR team as apparently 
asteroid but was later found to be a comet by other observers. Later on Smithsonian 
Observatory astronomers (IAU Cir. No. 7546, 2001) successfully photographed the comet 
with 1.2-m reflector of Mount Hopkins, Arizona and at that time it’s coma was 10 arcsec 
in size. 


On February 1, 2002 a comet C/2002 Cl was discovered visually by veteran comet 
hunter Kaoru Ikeya with a 25-cm aperature reflector and by Daqing Zhang with a home 
made 20-crn reflector. Ikeya has discovered four of the brightest comets of the past 
century, and this is his sixth comet. Zhang has been searching for comets for many years 
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, but this is the first comet to carry his name. Marsden and Nakano (IAU Cir. No. 
7843, 2002) have noted the similarity of the orbits of comets C/2002 Cl and C/1661 Cl 
and the numerical integration of the 2002 orbits backwards yields a previous perihelion 
date within a couple of years of 1659, making the link rather likely. During March and 
April 2002, comet C/2002 Cl was observable to northern hemisphere observers low in the 
western sky. It moved northward in the sky till 25 April. It has been faintly visible to 
the naked eye from dark sky sites in March and remained a binocular object till late May 
or June. During March and April the comet showed an obvious tail which faded rapidly 
in April. Photographs taken on March 18, 2002 using CCD camera at the Cassegrain 
focus of 104-cm telescope of State Observatory Nainital showed well developed tail and 
emission of discrete jets. The heliocentric distance on this date was 76 million km and 
the geocentric distance was 120 million km. The tail was extended upto about 0.8 million 
km and coma diameter was around 92 thousand km. 

Extensive obervations of every comet are needed to determine production rates, col¬ 
umn densities and surface brightnesses. Time resolved imaging of both the plasma tail 
and coma are needed to investigate mechanisms giving rise to the observed morphologies 
and associated changes. Therefore we took spectrophotometric observations of comets 
C/2000 WM1 (LINEAR) and C/2002 C2 (Ikeya-Zhang) during December 2001, and 
March, April 2002 respectively. 


2. Observations and Data reductions 

The spectrophotometric observational system consists of a HR-320 spectrograph, a CCD 
detector, a detector interface and a computer. The spectrograph gives a dispersion of 
2.4 A/pixel with a grating having 300 g/mm blazed at 5000 A. At the entrance slit of 
the spectrograph, a circular diaphram of 9 mm diameter corresponding to 2.33 arcmin as 
projected on the sky (p) and centered on the coma of the comet was used. The IK x IK 
CCD chip covers about 2500 A wide spectrum in a single exposure. 

At least three spectral frames of the each comet were obtained every night in the blue 
and red part of the spectrum. Along with the comet, standard stars were also observed 
to calibrate the flux of the comet spectra. Sufficient bias, twilight flats and sky frames 
were also taken. Data reduction is done using spectroscopic reduction software package of 
IRAF. We have used bright A-type spectrophotomeric standards to calibrate wavelength 
of the spectra and hence our wavelength correction may have uncertainly of ~ 10A. The 
flux and wavelength calibrated spectra (converted to frequency) are shown in Figures 
1, 2 and 3. We have an usable range in spectrum from ~ 3550 A to 6000 A for all 
dates. However on 15, 16, 18 and 29 April 2002, we have observed in red region also upto 
7000A and combined the spectra after wavelength and flux calibration. Comet C/2002 
Cl is observed on six nights of Mar and Apr 2002 and comet C/2000 WM1 was observed 
on three nights of Nov and Dec 2001 (Nov 20, Dec 1 and 6). However, it was found to 
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have no significant spectral features on Nov 20. Basic data at the time of observations 
for both the comets are given in Table 1. 


3. Emission bands, column densities and production rates 

The prominent features, as can be seen in Figures 1, 2 and 3 are CN(Av = 0) at 3883 A, 
C 2 (Av = +1,0,—1) at 4695, 5165 and 5538 A respectively. These features are weaker 
for C/2000 WM1. In order to measure fluxes in these emission bands, the continuum 
in spectrum was located by selecting wavelengths 3600 A, 4350 A, 4850A, 5400A and 
5800A which are free from emission bands. The area of strong emission bands were 
measured and converted into the total flux. Emission band fluxes relative to C 2 (5165 A) 
is given in Table 2. 

The number of molecules of each species, contained in a cylinder of radius defined 
by the diaphram used, and extending entirely through the coma was evaluated using the 
expression (Millis et. al. 1982) which is given as, 


log M(p) = log F(p) + 27.449 + 2 log(Ar) — log g (1) 

where F is the observed flux in cgs units, p is the projected radius of circular region in the 
sky at A, r and A are the heliocentric and geocentric distances of the comet respectively 
in AU and g the fluorescence efficiency (in cgs units) per molecule at 1 AU. We used the 
values of fluorescence efficiency for C 2 and C 3 from Sivaraman et al (1987). Because of 
the swing effect g varies significantly for CN with the comets heliocentric radial velocity. 
In order to determine radial velocity, the orbital elements for C/2000 WM1 and C/2002 
Cl were taken from British Astron. Assoc. Cir. No.780, 2001 and British Astron. Assoc. 
Cir. No. 782, 2001 respectively. During our observations the radial velocity varied from 
13.5 km s ^ 1 to 30.0 km s ^ 1 for C/2002 Cl and was ~ 28 km s _1 for the C/2000 WM1. 
The value of g was obtained for all the values of radial velocities from the Figure of 
Tatum and Gillespie (1977). Therefore the variation of g has been incorporated in the 
calculation of column densities of the molecule CN. The column densities thus obtained 
are listed in Table 3. 

The column densities thus calculated were converted into production rates (Q), as¬ 
suming a Haser model, through the relationship given by A’Hearn and Cowan (1975), 


M{p) = QV~ l p 


K 0 {y)dy + (l/x)(l - 1/p) + Ki(px) - K x {x) 


( 2 ) 


where V is the velocity of released species; x is the ratio between p and daughter molecule 
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scale lengths; /i is the ratio between daughter and parent molecules scale lengths; Kq and 
K\ are modified Bessel functions of the second kind of order 0 and 1. Following Delsemme 
(1982) we assumed V = 0.58/v^. The parent and daughter molecular scale lengths are 
taken from Cocharan (1985). The resulting production rates are given in Table 3. 

Production rates for dust have been estimated on the simplest possible model, spher¬ 
ically symmetric, uniform outflow ignoring variations in particle sizes, scattering angle, 
etc. The relationship derived by A’Hearn and Cowan (1975) is used to evaluate the 
production rate, Q , of dust in arbitrary units. 


Q = Lr 2 /p (3) 

where L is the luminosity of the comet at A = 4850 A. The dust production rates evaluated 
are listed in Table 3. 


4. Discussions 


The prominent emission features as seen in Figure 1, 2 and 3 are CN(Av = 0) at 3888 
A and C 2 (Ac = +1,0,—1) at 4690, 5160 and 5530 A. The strongest feature in the 
whole spectrum is due to C 2 (Ac = 0) at 5167 A. Figure 4 shows the extended region of 
the C/2002 Cl spectra around Na I emission at 5893 A, which has been detected also in 
the earlier observation (IAU Cir. No. 7914, 2002). However, Na I emission line was not 
taken into account for the purpose of calculating the parameters like production rates 
and column densities, as there was not sufficient signal at that wavelength. C/2002 Cl 
was nearest to the sun on March 23. Therefore the activity was maximum on that date. 
This is reflected in the continuum and in the production rates of molecules and dust 
which is maximum on March 23. As the heliocentric distance increased the activity was 
decreased. 

Figure 5 shows the production rates as a function of heliocentric distance for the comet 
C/2002 Cl and it varies approximately as r~ 2 for C 2 , CN and dust which is in agreements 
with the classical model of equilibrium vaporization of the nucleus. In the case of C/2000 
WM1 (LINEAR), the emission bands are weak although its total integrated magnitude is 
brighter. This may be due to its large heliocentric distance (at the time of observations) 
and the contribution of dust production is higher in comparison to others. The weak 
emission bands observed in C/2000 WM1 inspite of its bright integrated magnitude is 
probably because of the comets large heliocentric distance and because of much more 
dust content in comparison to other comets. 
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Table 3. Column densities ( M ) and production rates ( Q ) of the comets for the observed dates 
of Table 1. 


Comet 


log(M) 



log(Q) 



CN 

(Av = 0) 

c 2 

(Av = 1) 

c 2 

(Av = 0) 

o 2 

Av =-l) 

CN 

(Av = 0) 

c 2 

(Av = 0) 

dust 

4850A 

C/2000 WM1 

29.064 

28.347 

29.046 

28.429 

24.735 

24.927 

10.209 

(LINEAR) 

29.294 

27.946 

29.322 

29.042 

24.951 

25.192 

10.436 

C/2002 Cl 

31.165 

30.853 

31.108 

31.002 

26.706 

26.555 

11.693 

(Ikeya-Zhang) 

30.959 

28.889 

30.742 

30.636 

26.523 

26.514 

11.625 


30.841 

29.117 

30.655 

29.941 

26.404 

26.428 

11.265 


30.706 

29.825 

30.392 

30.509 

26.281 

26.176 

10.784 


30.629 

28.756 

30.300 

30.427 

26.202 

26.085 

10.848 


30.279 

28.841 

29.812 

30.139 

25.848 

25.591 

10.636 
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Figure 2. Absolute flux distribution of the head of the comet C/2002 Cl (Ikeya-Zhang) 



Figure 3. Absolute flux distribution of the head of the comet C/2002 Cl (Ikeya-Zhang) 
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Figure 4. The extended portion in the Na I emission region for the comet C/2002 Cl. 



Figure 5. The production rates of molecular species CN, C 2 and dust as a function of helio¬ 
centric distance for the comet C/2002 Cl. In the lower panel solid line represent CN and the 
dotted line represent C 2 






